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Abstract 

We present a self-contained formulation of spin-free nonrelativistic quantum mechanics that 
makes no use of wavefunctions or complex amplitudes of any kind. Quantum states are represented 
as ensembles of real-valued quantum trajectories, obtained by extremizing an action and satisfying 
energy conservation. The theory applies for arbitrary configuration spaces and system dimen- 
sionalities. Various beneficial ramifications — theoretical, computational, and interpretational — are 
discussed. 
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I. INTRODUCTION 



For nearly a century, quantum mechanics has presented philosophical and interpreta- 
tional conundrums that remain as controversial as ever. Far from disappearing into the 
realm of esoteric academic debate, recent experimental advances, e.g. in entanglement, de- 
coherence, and quantum computing, have brought such questions to the forefront of topical 
interest. The various competing viewpoints — Copenhagen,- Bohmian, 2-4 Many Worlds,- 
etc. — differ substantially in terms of their ontological interpretation of \I/ and its collapse, 
yet are all alike in their reliance upon a complex-valued wavef unction, and its propagation 
via the time- independent /dependent Schrodinger equation (TI/DSE). Evidently, a com- 
pletely self-contained, classical-like, and real-valued formulation of quantum theory without 
wavefunctions — if such a thing were possible — would necessarily present novel and poten- 
tially very important interpretational and computational ramifications. 

The aim of this paper is to show that for non-relativistic spin-free quantum mechanics, the 
above can be realized (A different approach is described in Ref. la) In particular, a quantum 
state can be represented as an ensemble of real-valued quantum trajectories, satisfying a 
self-contained partial differential equation (PDE). This was shown in a recent paper- by one 
of the authors (Poirier) for both the one-dimensional (ID) TISE and TDSE. It transpired 
that similar work had already been done for the TISE in ID,- the TDSE in 3D,- and in 
greater generality— In the current paper, we simplify, unify, and generalize these previous 
constructions, presenting quantum trajectory PDEs for arbitrary configuration spaces and 
system dimensionalities. The goal here is to present an alternative, standalone reformulation 
of quantum mechanics, that neither relies on the TDSE nor makes any mention of any 
external constructs such as and which, in addition, is likely to provide far-reaching benefits 
for numerical calculations, e.g. of accurate quantum scattering dynamics for chemically 
reactive molecular systems. 

The relevant PDEs can be derived from an action principle. We discuss symmetries of the 
action principle, associated conserved quantities, and other properties, such as a heretofore 
unexpected Hamiltonian structure in the case of ID time-independent quantum mechanics 
(TIQM), which also serves as the basis of an accurate many-D numerical scheme. We single 
out a specific choice of Lagrangian or gauge, and show why this choice may be regarded as 
physically preferred. Likewise, we single out a specific choice of trajectory-labeling coordi- 
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nate(s), in terms of which the resultant PDE exhibits no explicit coordinate dependences. 
We present simple analytical solutions, and discuss some initial numerical results that ap- 
pear very promising for molecular and chemical physics applications. The reformulation 
of quantum mechanics in terms of trajectory ensembles, in addition to shedding light on 
complex theoretical issues, evidently also provides important practical benefits. 

II. THE ID TIME-INDEPENDENT CASE 

As shown in Ref. 0, ID TIQM states can be represented uniquely with a single trajectory, 
x(t). The theory is one of a broad class of dynamical laws, for which the Lagrangian and 
energy are of the form 

L(x, x, x, . . .) = -rax 2 — V(x) — Q(x, x, . . .) , (1) 
E(x,x,x,...) = ^mi 2 + V(x) + Q(x, x, . . .) . (2) 

Equations (CEJ) and (J2J) are natural generalizations of the well-known classical forms. The 
quantum correction, Q, is similar to the potential, in that it appears with opposite signs in 
L and E, but is actually related to the TISE kinetic energy operator.- It has a universal 
"kinematic" form, i.e. no explicit x dependence. 

As in classical theory, the quantum trajectories, x(t), are obtained via extremization of 
the action, S — J Ldt. Since L is autonomous, the resultant x(t) solutions exhibit time- 
translation invariance and energy conservation for any choice of Q. However, for general 
Q, the conserved energy for Eq. ([I]) obtained via Noether's theorem does not take the form 
of Eq. (|2J). Requiring this equivalence imposes rather special conditions on Q, e.g., that Q 
must be invariant under time rescaling. The simplest, well-behaved, nontrivial Q giving the 
form of Eq. (j2j) is 

^/ "'\ h 2 ( x 5x 2 \ 

G(w) = i^-2iiJ- (3) 
Here fi is, in principle, an arbitrary positive constant. However, with the usual identifica- 
tion of h as Planck's constant, Eq. ([3]) is equivalent to the quantum potential of Bohmian 
mechanics^- for the ID Cartesian TISE. 7 This is surprising, given that the expression 
Eq. ([3]), being universal and kinematic, is determined entirely by the trajectory, x(t). In 
particular, no reference to the wavefunction, or the TISE itself, is used in this derivation. 
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Our quantum trajectories are nevertheless Bohmian trajectories, although our formulation 
and interpretation are not at all that of Bohmian mechanics, because no $ is involved. 

Action extremization applied to Eqs. ([T]) and ([3]) yields the following fourth-order au- 
tonomous ODE, describing ID TIQM quantum trajectories: 

dVix) h 2 fx 'x x , n x 3 \ n , AS 
mx + — ^ + — — - 8— + 10— = . 4 
ox 4m yx 4 x° x° J 

Two of the four constants of integration correspond to time-translation invariance and en- 
ergy, as in the classical case. The other two constants determine which particular TIQM 
state the trajectory is associated with. A one-to-one correspondence thus exists between 
trajectory solutions [of Eq. fl4|] and (scattering) TIQM states.- 

Any fourth-order ODE can be rewritten as a set of four coupled, first-order ODE's. 
Remarkably, Eq. (J3J) can be rewritten: 

s . dV(x) . 4rs 4 . p — s 8r 2 s 3 

x = i P = 5 i s = 72 i r = TT i (5) 

m ox mn m mn 



which are Hamilton's equations for a 2D system. 

OH . _ dH 

~o ' P n ' o ' " i 

op ox os or 



OH OH dH dH 

x = TT ' P = 5 r = ' s = — H - > ( 6 ) 



for the Hamiltonian 



tti \ s(2p-s) ... , 2r 2 s A 

H(x,p,r,s) = ^ J - + V(x) w . (7) 

2m mn 



In the above equations, (x, p) are the "classical" dimension phase space variables, and (r, s) 
correspond to an additional, "quantum" dimension, essentially describing quantum interfer- 
ence. Note that Eq. (J7|) reduces to the classical Hamiltonian when r = and s = p. 
In terms of the time derivatives of x(t), 

h 2 x h 2 fx 2x 2 \ 

S = mX]r = ^ ]P = mX + ^{x^-^)- (8) 

Substitution of Eq. OH]) into Eq. (JJJ) then reveals H to be the conserved Noether energy of 
Eq. (j2J). For the free particle case [dV/dx = 0], p is a second conserved quantity, in involution 
with H (and also derivable from Noether's theorem). The importance of p is difficult to 
overstate; it represents the "particle momentum," analogous to the well-known "particle 
energy," E = iJ.-^ Yet remarkably, p has barely been considered^ in the previous literature, 
which generally regards s = mx as particle momentum. Note that p is conserved for all free 
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particle TIQM states, including those exhibiting interference, whereas s is conserved only 
for plane wave states — i.e., the classical special case for which p = s. 

The above Hamiltonian approach is proving extremely useful in numerical calculations 
of quantum reactive scattering phenomena-i For any ID TIQM application, exact quantum 
results are obtained simply by propagating a single x(t) trajectory using the classical Hamil- 
ton's equations of Eq. (Q, until pit) flattens asymptotically to the desired level of numerical 
accuracy. The final p value then provides a direct measure of the quantum reaction probabil- 
ity. Note that Eq. ([6]) is amenable to efficient, symplectic numerical integrators as are used 
in classical simulations, and the conservation of E can be monitored as an on-the-fly measure 
of computed accuracy. This approach is extremely robust, accurate, and efficient, leading 
to 15 digits of accuracy, even in the extremely deep tunneling regime (where absorbing po- 
tentials render conventional exact quantum scattering methods intractible). For many-D 
TIQM applications, classical-like sampling over quantum trajectory initial conditions leads 
to an approximate simulation scheme that has also proven to be remarkably accurate (i.e., 
to two or three digits). 11 



III. THE ID TIME-DEPENDENT CASE 



For the case of ID time-dependent quantum mechanics (TDQM), any self-contained for- 
mulation must involve a PDE, rather than an ODE. It is no longer possible to exactly 
represent a quantum state as a single trajectory, x(t), but rather as a one-parameter ensem- 
ble of trajectories, x(C,t), where the real- valued, space- like coordinate C labels individual 
trajectories. The equation of motion should be a PDE involving C and t derivatives, prefer- 
ably derived from a field-theoretic action principle. 

In Refs. , and llOl. C was chosen as the initial trajectory value [x(C, 0)=Xq=C]. The 



resultant PDE is complicated, exhibits explicit xq dependence through the initial probability 
density, p (x )=p(x , 0), and bears little resemblance to Eq. In addition, Q is expressed 
in terms of C rather than t derivatives of x. The PDE can be simplified by a better 
choice of the trajectory parameter, C, which in general can be taken to be any monotonic 
function of xq (regardless of the initial wavefunction) . A crucial idea of the current paper 
is that C should be chosen so as to uniformize the probability density. In particular, since 
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pc{C) dC = p(x, t) dx, if we choose 



C 



j-t) 



po(x' ) dx' 



(9) 



then C takes values from to 1 (for normalized wavepackets), and pc{C) = 1. 

Working with Eq. (or any uniformizing choice of C), and writing x' = dx/dC, x 
dx/dt etc., the PDE of Ref. l7| simplifies very substantially to 



.. dV(x) h 2 fx"" 
dx Am \ x' 4 
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(10) 



Equation ( ITUj) is the perturbed Newton equation for the ID TDQM case; it has no explicit 
coordinate dependences. It also bears an extremely close resemblance to Eq. (JIJ), obtained by 
replacing C derivatives with t derivatives in the last term on the left hand side (representing 
the quantum force). More formally, Eq. (j3J) is obtained on looking for travelling wave 
solutions of Eq. ({TO]) , i.e. solutions of the form x(C, t) = x(t — AC) where A is a constant. In 
the ID TIQM context, t thus serves as an effective uniformizing coordinate. Equation ( flOl) 
also admits t-independent solutions, which correspond to the bound (ffuxless) ID TIQM 
quantum states (the previous section concerns only the scattering states). 

As in the ID TIQM case, Eq. ( II 01) is a variational equation, obtained by extremizing the 
action 



dCdt 



-mx 2 — V(x) | ~~o 
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(11) 



c.f. Eqs. ([T]) and ([3]). This action is invariant under translations of both coordinates t and 
C . By Noether's theorem this gives rise to two conservation laws, which are easily found to 
be, respectively, 
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The first corresponds to conservation of energy. In the free particle case, there is also a 
momentum conservation law, arising from ^-translation symmetry: 



d r . h 2 d 
M [mx] + Am-dC 
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Interpreted as hydrodynamical balance equations of the general form OA/ dt+dB / dC = 0, 
the first square bracket in Eq. ffl2|) [Eq. ffT4]) ] represents the energy [momentum] density, and 
the second term the corresponding flux. This designation is only determined up to addition 
of a C-derivative to the density A and subtraction of the corresponding t-derivative from 
the flux B. Consequently, the energy density need not conform to the standard TDQM 
"field" form,- but instead may be chosen to be the particle energy, T+V+Q, as in Eq. (fT2l) . 
This choice is appropriate for the Lagrangian T— V— Q [cf. Eqs. ([I]) and (j2J)], and has 
the great advantage of being conserved along individual trajectories in the TIQM limiting 
case (in general, only the total ensemble energy is conserved). However, in the momentum 
conservation law Eq. fTH|) . it does not seem to be possible to use a density that reduces to 
the particle momentum p in Eq. ([S]) • 

Also of note is the balance equation, Eq. (fTSl . For any autonomous Euler-Lagrange PDE, 
the Lagrangian density L is determined only up to the addition of a divergence (i.e., the sum 
of C- and t-derivatives). We have chosen forms of the action [Eq. ( Till ] and the conservation 
laws such that L appears as (minus) the flux in Eq. ( 1131) . Via gauge transformations, it is 
possible to eliminate the third-order derivative from the Lagrangian and energy densities, 
but not from the flux of Eq. (fT3l) . The choice of L we have made has the advantage that 
the trajectory action, S(C,t) = J * L(C,t')dt' ', expressed in units of %, can be identified with 
the change in phase of 

We now consider Gaussian wavepacket evolution under the free particle (V = 0) and 
harmonic oscillator (V = \moj 2 x 2 ) potentials. The respective x(C, t) solutions are 

x + + a erfinv(2C - 1)^1 + and (15) 

x cosu(t - t ) + posi "^' fo) + a erfinv(2C - 1) x ^cos 2 u(t - t ) + * gg^ to) , (16) 

where x ,p ,t ,a are real wavefunction parameters, and < C < 1. Note these solutions 
diverge as C — > or 1. In general, x(C, t) must diverge at the C endpoints; modulo this 
requirement, any solution of Eq. f flUj) (or its arbitrary-C generalization) can be used to 
reconstruct a normalized solution \I/(x, t) of the ID TDSE. 

We have successfully numerically integrated Eq. ( JTOj) for an Eckart potential and initial 
Gaussian wavepacket, using the Stormer-Verlet and other methods. Initial results compare 
favorably with those obtained using standard \l/-based methods, but further improvements 
in efficiency and accuracy are planned. In any case, these calculations represent a milestone 
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achievement, as the first successful synthetic Bohmian quantum trajectory calculations ever 
achieved for a system with substantial reflection interference — a much-sought goal eluding 
chemical dynamics researchers for over a decade, and a major hurdle preventing exact quan- 
tum wavepacket calculations for large molecular systems with few reaction pathways.- 1 ^ 



IV. THE MANY-D TIME-DEPENDENT CASE 



The ID analysis generalizes to many-D. The single variable x is replaced with the n- 
dimensional configuration space vector x, with C likewise replaced with C, so that x(C, t) 
represents an n-parameter family of trajectories. One option^ is to take C = x , though 
we wish to consider more general choices for which C and xo are related via any invertible 
coordinate transformation. 

In analogy with the ID case, it can be shown that for an arbitrary choice of the 
parametrization C, 

where J is the Jacobi matrix, = dx % jdC 3 . As in the ID case, we mostly work with 
a uniformizing C for which pc(C) = 1. The resulting perturbed Newton equation, i.e. 
the many-D generalization of Eq. ( flOl) . can be written in various different forms, the most 
compact being 

dx l 4mdC m \ 1 3 dC k dC l ) v ' 

Here K = J~ l denotes the inverse Jacobi matrix. The Einstein summation convention 
is used, albeit with some mismatched indices as we are currently assuming the Euclidean 
metric on x space. 

Equation ( TTHl) is the variational PDE for the action 

1 



al n Calt 



m x • x — V(x.) — Q 



where (19) 



h 2 ( k d 2 K) \dK\dK k \ 
Q ~ ~^m~ [ Kj dC l dC k + 2 dC l dC k ) ' (20) 

Equation f lT9"|) is the many-D generalization of Eq. ffTT]) . which preserves the L = T — V — Q 
and \1/ phase properties. Like Eq. (11 II) . Eq. (1191) has associated laws of conservation of 
energy and momentum. Invariance under C translation is replaced by invariance under the 
infinite-dimensional group of volume preserving diffeomorphisms. 
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For a general (not necessarily uniformizing) choice of C, and allowing a non-Euclidean 
metric g%j{x) on x-space [note in this case a factor of y/g should be inserted in the denomina- 
tor of the RHS of Eq. (JT71)]. the action, in a simpler gauge involving only second derivatives, 
becomes 



J J pc(C) d n Cdt 



-mg ij (x) x l xi - V(x.) 



8m 
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l dC k \ det J ) dC l \ det J 



(21) 



cf. Eq. (3.5) in Ref. |_10|. Note that a many-D solution of the TDSE gives rise to a solution 
of Eq. (118p . but the converse holds only if the initial x field is a gradient; this property is 
then preserved by the evolution Eq. ( TT8"]) . 10 



V. CONCLUDING REMARKS 



We have developed a self-contained, trajectory-based formulation of spin-free nonrela- 
tivistic TDQM, achieving all goals as outlined in the Introduction. Further developments 
are underway. Theoretical progress will require a correct treatment of spin, relativity, parti- 
cle indistinguishability, and second quantization — with a promising start having been made 
by Holland and others.-"^ The invariance of Eq. ( I2TT) under volume-preserving diffeomor- 
phisms suggests a connection with gravity, though the physical significance of C is not yet 
clear. 

Numerically, the prospect of stable, synthetic quantum trajectory calculations for many-D 
molecular applications will be fully explored, as the benefits here could prove profound.-^ 1 ^ 
Our formalism offers flexibility for restricting action extremization to trajectory ensembles 
of a desired form (e.g., reduced dimensions), thereby providing useful variational approx- 
imations. Our exact TDQM equations are PDEs, not single-trajectory ODEs — the entire 
ensemble must be determined at once. But they provide the great advantage of making 
no reference to any external fields, such as densities or wavef unctions. Alternatively, the 
many-D Hamiltonian ODE approach is approximate, but evidently quite accurate.— 

Regarding interpretation, we draw no definitive conclusions here. However, it is clearly of 
great significance that the form of Q can be expressed in terms of x and its C derivatives — 
implying the key idea that the interaction of nearby trajectories, rather than particles, is the 
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source of all empirically observed quantum phenomena (suggesting a kind of "many worlds" 
theory, albeit one very different from Ref. (5j). As such, it is locality in configuration space, 
rather than in the usual position space per se, that is relevant. In effect, we have a hidden 
variable theory that is local in configuration space, but nonlocal in position space — though 
the latter is hardly "spooky" in the present nonrelativistic context [even classical theory 
is nonlocal in this sense, depending on U(x)]. Many ramifications are anticipated for a 
wavefunction-free interpretation of measurement, entanglement, etc. One wonders whether 
Bohm would have abandoned pilot waves, had he known such a formulation was possible — 
or, for that matter, whether the notion of quantum trajectories might have actually appealed 
to Einstein. 
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